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Abstract 

We investigate in some detail consequences of the effective colour-dielectric formulation of 
lattice gauge theory using the light-cone Hamiltonian formalism with a transverse lattice. As a 
quantitative test of this approach, we have performed extensive analytic and numerical calcu- 
lations for 2 -I- 1-dimensional pure gauge theory in the large N limit. Because of Eguchi-Kawai 
reduction, one effectively studies a 1 -I- 1-dimensional gauge theory coupled to matter in the 
adjoint representation. We study the structure of coupling constant space for our effective po- 
tential by comparing with results available from conventional Euclidean lattice Monte Carlo 
simulations of this system. In particular, we calculate and measure the scaling behaviour of the 
entire low-lying glueball spectrum, glueball wavefunctions, string tension, asymptotic density of 
states, and deconfining temperature. We employ a new hybrid DLCQ/wavefunction basis in our 
calculations of the light-cone Hamiltonian matrix elements, along with extrapolation in Tamm- 
Dancoff truncation, significantly reducing numerical errors. Finally we discuss, in light of our 
results, what further measurements and calculations could be made in order to systematically 
remove lattice spacing dependence from our effective potential a priori. 
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1 Introduction 



The least understood scale of the strong interactions is the intermediate one, between short distances 
described by asymptotically free QCD and the relatively long range phenomena of nuclear forces. 
A successful description of this intermediate region should not only include 'global' properties of 
hadrons, such as masses and decay constants, but also the wealth of accumulated experimental 
data on hadronic sub-structure. In addition one might hope to shed some light on the interface 
between gauge theory and nuclear physics. There are also many reasons why knowledge of hadron 
wavefunctions in a general Lorentz frame, particularly at the amplitude level, will be necessary 
for future progress in all aspects of particle physics. At present, no single approach to the strong 
interactions seems capable of this. In this paper we pursue an idea which, in principle, can address the 
above problems. As we will demonstrate, the colour-dielectric concept is both a natural description 
and, when combined with the light-cone Hamiltonian framework, a practical description of the strong 
interactions at intermediate scales. A number of authors have proposed the colour dielectric idea in 
various forms [Q-Q^but, while the motivation was often elegant, results were usually of a qualitative 
nature. Our work is strongly foreshadowed by the work of Bardeen, Pearson, and Rabinovici [Q in 
which the 'dielectric' formulation of a light-cone transverse lattice gauge theory was first proposed 
(although the authors did not discuss it in quite those terms). 

We characterise a dielectric formulation as one in which gluon fields, or rather the SU{N) group 
elements they generate, are replaced by collective variables which represent an average over the 
fluctuations on short distance scales. These dielectric variables carry colour and form an effective 
gauge field theory with classical action minimised at zero field, meaning that colour fiux is expelled 
from the vacuum at the classical level. The price one pays for starting with a simple vacuum 
structure, which may arise only for a rather low momentum cut-off on the effective theory, is that 
the effective action will be largely unknown and must be investigated per se. 'Conservation of 
difficulty' appears to have dashed one's hopes. However, there are at least two reasons why shifting 
the complication from the vacuum to the action is desirable. Firstly, although most approaches to 
field theory place great emphasis on the calculation of vacuum energy and field condensates, neither 
of these are measured in experiment. Secondly, if one hopes to employ light-cone quantisation, which 
seems the only practical way to calculate the hadron's wavefunction in a general Lorentz frame, then 
dealing with explicit vacuum structure can be rather awkward. 

We choose to employ a lattice dielectric gauge theory in which the fields are by complex 
matrices M, in the adjoint representation of the gauge group, attached to links of the lattice 

■^A review of subsequent investigations and more extensive references can be found in Ref. |j7|. 



2 



Because the lattice spacing in this theory should be rather large if the colour dielectric idea is to be 
valid, it is essentially impossible to obtain the effective action by explicitly integrating out gluons 
(blocking from the continuum). More realistic, but still essentially impossible, would be to introduce 
the dielectric fields by blocking the unitary link fields U E SU {N) of Wilson's lattice gauge theory 
That leaves a 'trial-and-error' approach where one tunes the coefficients of the effective action, 
or some truncation of it, until solutions exhibit the correct symmetries, such as Lorentz invariance. 
This more-or-less rules out a Euclidean path integral formulation, since one would have to do the 
integral (numerically) for each set of coefficients. A Hamiltonian approach is a viable alternative 
since most of the work goes into computing matrix elements of the Hamiltonian. This only needs 
to be done once for each operator appearing in the effective action. It is only the diagonalisation of 
the Hamiltonian, which can be performed much more rapidly, that needs to be repeated for different 
coefficients of the operators. In the light-cone Hamiltonian formulation the continuum limit is 
taken in the temporal and in one spatial direction, leaving a transverse lattice. By using a light- 
cone Hamiltonian formulation the quantum vacuum remains as simple (modulo zero modes) as the 
classical vacuum. 

Next we must choose a specific problem in order to test the practicality of this dielectric formu- 
lation. In this paper we postpone the introduction fermions (and hence quark structure functions) 
and study pure gauge theory. We look at glueballs in three spacetime dimensions at large N. 
The glueball problem is one of the most difficult one could imagine since it is ultra-relativistic and 
strongly interacting. However in four dimensions there is a dearth of accurate, detailed experimental 
or theoretical results which could be used as a benchmark. While there is an undoubted urge to 
contribute to the presently controversial debate on glueballs in nature, as a first calculation a more 
circumspect choice is appropriate. Fortunately there is a wealth of accurate information now avail- 
able on glueballs for pure gauge theory in three dimensions from the recent Euclidean Lattice Monte 
Carlo (ELMC) simulations of Teper ||. That work has shown that pure non-Abelian gauge theory 
behaves much the same way in three as in four dimensions: a discrete set of massive boundstates 
are generated by a linearly confining string-like force. Moreover, Teper has performed calculations 
for 2, 3, and 4 colours, allowing an extrapolation to the large N limit. The relevant expansion 
parameter is and observables measured at ~ are extremely close to those measured at 

— 1/9. The large- iV limit is convenient, though not essential, for our Light-Cone Transverse 
Lattice (LCTL) formulation, which becomes simpler in the limit of many colours. It enables us to 
implement Eguchi-Kawai reduction [ pO[ to a 1 -I- 1-dimensional theory which does not need to be 
quenched or twisted so long as the vacuum is the simple dielectric one. It also allows us to make use 
of the extensive intuition developed for large-A^ theories in relation to string models. Indeed, our 
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formulation offers the rare possibility of describing the parton, constituent, and string behaviour of 
hadrons in one framework. The relationship between these pictures, each very different but equally 
successful, remains one of the outstanding enigmas of QCD. 

We should emphasize that this is the first time that the LCTL formulation has been used in a 
serious quantitative study, and the expectations of the reader should be tempered accordingly. It 
was necessary to overcome a number of technical problems whose solutions probably have wider 
applicability. Our initial goal has been to explore the coupling constant space of the dielectric 
effective potential, which was largely ignored by Bardeen et alii in their pioneering study. We 
will fix our effective potential by comparison with the the ELMC spectrum. This gives rise to 
predictions for other observables and it is non-trivial that a consistent picture emerges. In the light 
of our results, we will also discuss ways of measuring further observables to test Lorentz covariance 
and so fix the effective potential independently, as well as the sine qua non of a systematically 
improvable analytic renormalisation group analysis. 

The outline of the paper is as follows. In the following section we describe in general terms 
the dielectric formulation on a Euclidean lattice, with particular reference to the large N limit. 
The associated light-cone Hamiltonian limit of the transverse lattice is constructed in Section ||. In 
Section 1^ we obtain analytic and numerical solutions for glueballs, compare with Teper's data 
and identify a scaling trajectory in coupling constant space. Section |] addresses the string tension 
via winding modes. We discuss the asymptotic density of glueball states and associated finite 
temperature properties in Section |^ and compare with results from string theory. The conclusions 
in Section ^ contain discussion of some general points of principle, as well as the shortcomings of 
LCTL as used in this paper and their possible resolution. Appendix ^ contains some details of 
the proof of Eguchi-Kawai reduction in light-cone Hamiltonian formalism. In order to obtain the 
required numerical accuracy, we had to devise new methods of computation and extrapolation for 
light-cone Hamiltonians. These important techniques are collected in Appendices B and C since 
they are relevant to many parts of the paper. 

Some of our preliminary results have been presented in Refs. 

2 Dielectric gauge theory on a Euclidean lattice 

In this section we motivate the colour dielectric formulation of gauge theory in the context of a 
Euclidean lattice ||^. For intermediate scales, we stress that we have not derived it in detail from 
QCD at short distances, but view it as a starting point, similar to a sigma model for mesons. 
However, unlike a sigma model, or even a quark effective theory, there will be explicit colour-field 
degrees of freedom. The best we can do at this stage is to discuss the general form expected of the 
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effective potential. The details for the transverse lattice will be determined later by actually solving 
the theory as a function of the associated coupling constants. 

2.1 Effective action 

Let us consider Wilson's formulation |^ of lattice gauge theory. The principle observation is that, 
for large lattice spacing, the unitary link variables U have large fluctuations away from the identity 
and it may be just as well to replace them by some other set of variables which takes advantage of 
this randomness. Nielsen and Patkos j|] chose a 'block' transformation from the product of parallel 
transporters U along some set of lattice paths Cat weighted by p{Cab) between two points a and b of 
a sub-lattice 

M{a,h) = p{Cab)\{U . (1) 

{Cab] Cab 

The blocked matrix link variable M G GL{N,C), since it is a linear combination of fairly random 
unitary variables. This blocking scheme has been investigated by operator matching using Monte 
Carlo and Schwinger-Dyson techniques j7j for SU{2), with conclusions for the general form of the 
effective potential similar to those described below for large N. Although we will not use the 
reasoning (|l|) directly, it motivates us to search for an effective lattice gauge theory in terms of such 
non-compact link variables. 

One knows that at large lattice spacing a Wilson's action gives a fair description of non-Abelian 
physics. One assumes that as a — s- the lattice errors gradually reduce to nil. In order to reduce 
the lattice errors at a finite a one might imagine 'smearing' the link variables U. Therefore consider 
introducing both a unitary matrix U{1) and a complex matrix M{1) on each link I of the spacetime 
lattice, with partition function 

Z = J Y[VUil)VM{l) exp (^^E Ti-pM - y E ^^{(^(0 " t/(/))(Mt(/) - U^l))}^ (2) 

where Tr p [M] denotes the trace of the product of link variables M around an elementary plaquette 
P and gL is the dimensionless lattice gauge coupling constant. V denotes the usual Haar measure, 
which is simply 

N 

VM ^ W dM,j dM*. (3) 

for the complex matrices. In this case we have introduced a simple Gaussian smearing function 
which peaks the measure of the complex matrices on unitary matrices. In the limit A ^ 0, one 
recovers Wilson's theory. It is natural therefore to associate small A with small a. For small A one 
can derive a series of corrections to the Wilson action Trp[J7] as an expansion in A. Define the 
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complex matrix 



C(0 = ^/- [M{l)-Uil)] 



(4) 



If Ip, 2p, 3p, and 4p are the four links that form plaquette P, we have 
Z 



(5) 



where 



-Ay[c/] 



exp , „ , 

\9l / 



J llVUil) exp Trp[[/] j + ^Ay[C/]j 

n2?C(0exp (-^Tr{Ct(OC(0} 

+ E Tr {C(lp)C/(2p)C/(3p)[/(4p) + C/(lp)C(2p)[/(3p)(7(4p()6) 

51 P 

+C/(lp)C/(2p)C(3p)C/(4p) + t/(lp)C/(2p)C/(3p)C(4p)} + 0(A)) 

and we have dropped an irrelevant normalisation constant. Performing the Gaussian integrals over 
C{1) one generates an effective action for [/-variables consisting of arbitrary Wilson loops: 



Ttp[U] + AV[U] = ^4 
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(7) 

(Only the terms in one plane are shown.) The leading correction is the well known rectangular 
C/-loop which, with the correct coefficient, will cancel O(a^) errors of the Wilson action. One 
might imagine that by introducing a sufficiently complicated smearing function one could produce 
systematically the terms needed to correct for higher lattice spacing errors, id est a perfect action 
which tracks an infrared stable RG trajectory ||T^. With the simple Gaussian smearing, the first 
correction to Wilson's action (^ in fact has the wrong sign for improvement. However we are free 
to make the effective potential more sophisticated beyond the Gaussian approximation by explicitly 
adding zero-area Wilson Af-loops (which become trivial in the limit A ^ 0) to the action (^; 

, 2 



J2 T^{M\l)M(l)M\l)M(l)] , ^(Tr{Aft(Z)M(0})^ 

et cetera] 



Tt:{M{1)M{1')M\1')M\1)] , 

{IM) 

{I, V) denotes a sum over neighboring links. Note that gauge invariance 

M{1) -> Vl^ iMil)V,j+i , V e SUiN) 



(8) 



(9) 



provides a strong constraint on the allowed terms we can add at each order in the field M . Contri- 
butions to the effective potential of the form 



^ TT{M{l)M{l + a)M\l + a)Ml{l)} , 

A 



(10) 
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where > and ^ + a is the Hnk consecutive to / (parallel to I and sharing one site) , will be necessary 
in the small lattice spacing limit A ^ in order to have any chance of producing the correct sign 
for the coefficient of the leading rectangular CZ-plaquette correction: the coefficient in Eqn. (|^) now 
becomes A(l — 2v)/ g\N . Such terms have generally been neglected in previous studies of dielectric 
lattice gauge theory. 

We may instead integrate out the [/-variables link-by-link using the large- saddle-point result 
of Ref. il3|. Thus, for some link L 



j VU exp (^-^Tr{(M- [/) {M^ - U^)}^ = 



(N N 
-N^x - (AA, - 2VA7T^) - i ^ log (VA7T¥ + V^vjT^) - (! + x) 1 (11) 

where A, are the eigenvalues of MM^/A^ and 



1 ^ 1 

x = Q if — y^<2 



1=1 



ly—^ = 2 if ly^>2 (12) 

Note that A and A^ are 0(1) while the effective action is 0{N'^) and g\N is, as usual, held finite 
in the large A^ limit. The asymptotics of the effective M-potential ( pT| ) versus TrjA/^A/} /N are 
easily found: for large eigenvalues A^ we are in the x — Q regime and so ( pT| ) becomes 

exp f-y TrjikftAf} j , Tr { Af^Af} ^ oo ; (13) 



for small A^ we are in the x ^ regime and may expand in powers of A, 

1 1 

4 ^ iVA2 

to obtain 



^ = 7-l^Tr{AftM}+0(MVA^) (14) 



exp(^(^^ + ^^ ATr{A//tA/}+0(AfVA'')^ , ^ Tr {A//t Af} -> . (15) 

As A is reduced through 1 the (classical) minimum shifts from M = to TrjAftAf} > 0. We 
should not take too seriously the details of the potential derived from ( pT| ) beyond the Gaussian 
approximation, since the additional terms (|8|) must then be taken into account. 

We see the qualitative picture anticipated for the dielectric formulation emerging. At large 
enough lattice spacing (large enough A), the effective [/-potential (|^), which contains arbitrary 
Wilson loops with roughly equal weight, would be re-organised by the dielectric variables into a 
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small-field expansion about M = 0. In practice one will have to truncate the small-field expansion 
and the form of the effective M-potential should be fixed by scaling analysis. It is difficult to do 
this analytically because, apart from the mass of the M-field, there is no identifiably small coupling 
at hand; A > 1 is not small, so we cannot perform improvement near the fixed point, as is now 



popular in the ELMC approach |14 . In fact, the ultimate justification for expecting a truncation to 
work comes from light-cone quantisation, as we argue in the next section. (We will use an effective 
potential truncated to fourth order in the fields (|35|).) We anticipate that there should be an optimal 
'window' for use of such a truncated potential. A heavy M-field, when the minimum at M = is 
steep, presumably implies that one is deep in the dielectric regime, where the lattice spacing must 
be very large (analogous to the very strong coupling regime of Wilson's theory); in this case any 
truncation of higher order terms in may result in severe finite lattice spacing errors. On the 

other hand our Hamiltonian analysis will only be valid for quantisation about the M = classical 
solution, so the mass squared of the A/-field must be positive. This suggests we should take it to be 
positive and small. Our results will bear this out. 

Finally we note that one would arrive at basically the same conclusions for an Abelian theory. 
However it is well-known that there is a phase transition between strong and weak coupling in this 
case [||, The dielectric formulation in the regime where M has positive mass squared would 
correspond to strong coupling, and is thus not useful for studying weakly coupled QED. 

3 Light-cone Hamiltonian formulation 

3.1 Transverse lattice 

A Hamiltonian limit of the dielectric lattice theory requires us to take the lattice spacing to zero in 
at least one direction. We now specialise to 2 4- 1 dimensions, since it is for this case that we will 
perform explicit calculations. Since we also wish to make use of the fact that two components of 
the vector potential are non-dynamical, we take the continuum limit in the x'^ and directions, 
leaving the transverse direction x^ discrete. This is means that only dynamical degrees of freedom 
will be involved in the effective potential V^[Af], a simplification not shared by the Euclidean lattice 
or equal-time Hamiltonian quantisation of this problem. We now use different smearing parameters 
A and A for the transverse x^ and longitudinal {x'^,x'^) directions, respectively; likewise we define 
transverse and longitudinal lattice spacings a and a. For the longitudinal directions we take the 
limit A, a ^ so the measure becomes infinitely peaked on unitary matrices; this ensures that 

lim M{1) = U{1) = 1 + ihAa + ■■■ (16) 
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for a link I in the direction a = or 2. Applying Eqn. ( p^ ) to Trp[A/] one derives a transverse 
lattice action whose form was first suggested in Ref. ^ 

A = J da;0da;2^(^Tr{i^„M,ip"A40^}-4^Tr{F„^f^"/'} 



V,i[M]\ (17) 

J — \ -^y- J 

where 

DaM^i = (da + iAa{x^)) M^i - iM,j.iAa{x^ + a) (18) 

We have introduced the dimensionful continuum gauge coupling = ag\/d? in 2 + 1 dimensions 
and rescaled M MgL- Afj.i lies on the link between and x^ + a while Aa{x^) is associated 
with the site x^ . V^i [M] is a purely transverse gauge invariant effective potential, which is assumed 
analytic about its minimum at M = 0. For 2 + 1 dimensions this means it contains only (products 
of) zero-area Wilson loops such as those displayed in (|[) with one point pinned at x^. 

A feature of light-cone formalism, which is not present in the equal-time approach, is that for 
free fields of mass squared fi^ > the momentum fc+ ~ [k^ + fc^)/V2 is positive definite for finite 
energy since 

--^^ 

This, together with momentum conservation, implies that even the quantum vacuum is free of 
particles and that a large number of particles necessarily carry a large light-cone energy. It is for 
this reason that one believes a calculation of the low mass spectrum using the light-cone Hamiltonian 
is feasible when effects due to large numbers of particles have been truncated. Higher powers of the 
field M in the effective potential couple precisely to those components of the wavefunction with 
a large number of A/-quanta. Of course, in the LCTL construction, the Af-quanta are also link 
variables, so energy considerations are intimately tied in with questions of spacetime symmetry 
restoration. Note that the above arguments break down if applied to the gluon fields Ai directly 
since they are massless and can be produced copiously for = 0. On the other hand, if one gives 
the gluons Ai a fake mass this breaks gauge invariance leading to copious production of non-local 
terms in the effective potential. 

3.2 Gauge symmetry 

Having formulated a transverse lattice action ( |l7|) we now introduce light-cone co-ordinates x^ = 
Xzf: = (x° ± a;^)/\/2 and quantise by treating x+ as canonical time, following Ref. ||l[. The theory 
has a conserved current 

J+(x) = I [a//,! {D.M^i)^ ~ {D_M,i)mI + Ml_^D^M,i_a - {D^M^i^a)^ Af,i_J (20) 



Note that this definition of g differs from the standard one by a factor of 2. 
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at each transverse lattice site which plays a special role. If we pick the light-cone gaugeg A_ = 
the non-propagating field satisfies a simple constraint equation at each transverse site 

{d^f A+{x') = . (21) 

a 

Solving this constraint leaves an action in terms of the dynamical fields M^i 

J dx+dx- J2 Tr l^aM.ia^M], + |^ J+ j^^i | " [M] , (22) 

the pole in the propagator given by 1/(9-)^ is regulated by a principle value type prescription 
[Q, Q . It is straightforward to derive light-cone momentum and energy 

P+ = 2y"dT~^ Tr|a_A4ic)_Mti| (23) 
P- = j dx-Y^VAM]-^^jAj+^J+\ (24) 

There remains residual gauge symmetry under ^"-independent transformations at each site 

M^i , V G SUiN) (25) 

Evidently the light-cone energy (e3) will not be finite unless the associated charge vanishes at each 



dx'j\=0 (26) 



This requirement also follows from the zero momentum limit of the constraint equation (|21| ) and will 
force one to consider only singlet states under the transformations ( p5| ) . This condition is analogous 
to the colour neutrality imposed by Kogut and Susskind [T^ as a consequence of Gauss' law, but 
here it follows from dynamical considerations alone. 

3.3 The Hilbert space 

Expanding in harmonic modes at x~^ — about M = 0, it is convenient to work in longitudinal 
momentum space and transverse position space 

M,.{x-) = -^j^ ^[a.i{k,x')e-'''^' +{a+iik,x')y e'''^') (27) 

where the modes satisfy equal-x+ commutators 

a\,i.j{k,x'^), {ap^ki{k,x'^))'' 
ax,ij{k,x^),apM{k,x^) 



= SikSji6xpS^i,j;i6{k - k) (28) 
= (29) 



*We set the associated zero mode J dx A- to zero as a dynamical approximation. 
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with colour indices i,j e {1, . . . ,N}, orientation indices A, p G { — 1,+!}; {a\.ij) = {'^x)ji ^-nd for 
clarity we omit the + superscript on k. The Fock space operator a'^i{k, x^) creates a mode on the 
link between sites and + a with longitudinal momentum k and orientation ±1. In the Fock 
space, only combinations singlet under (^5|) are annihilated by the charge (|2^). This gives a Hilbert 
space at fixed x^ formed from all possible closed Wilson loops of link modes a± on the transverse 
lattice. Thus, a typical p-link loop will be something like 

Tr {al^{ki,x^)al,{k2,x^)al,{k3,x^ ~ a) ■■■ al^{kp,x^ - a)} |0) (30) 

The number of +l's equals the number of — I's for a closed loop on the transverse lattice, while the 
longitudinal momenta km are unconstrained except for X)m=i ^™ ^ ^^'^ ^™ — ^- abuse 
of language, we will often refer to p as the number of particles — the link variables are reminiscent 
of constituents — even though the relevant field variables are non-locally defined. At large N we 
need only study the dynamics of single connected Wilson loops in the Hilbert space — there are no 
terms with more than one trace — since the loop-loop coupling constant is of order These 
loops may be thought of as 'bare' glueballs, and the problem is to find the linear combinations 
that are on mass shell. Neglecting k^ = 0, which is consistent with expanding about the M = 
solution of the dielectric regime, the Fock vacuum is an eigenstate of the full light-cone Hamiltonian 
p- 10) = P+ |0) = 0. 

3.4 Eguchi-Kawai reduction 

A further simplification occurring at large N is Eguchi-Kawai reduction ||ic| ]. This means that in 
the Lorcntz frame = the theory is isomorphic to one compactified on a one-link transverse 
lattice with periodic boundary conditions, where P" acts on a basis of zero winding number Wilson 
loops ||l7|| , id est we can simply drop the argument x^ or I from M in all of the previous expressions. 
Effectively one is now dealing with a 1 + 1-dimensional gauge theory coupled to a complex scalar 
field in the adjoint representation (with self-interactions). Mathematically isomorphic "coUinear 
QCD models" have recently been studied ||2^, the equivalence following after identifying winding 
number (orientation) here with helicity there. 

The Eguchi-Kawai reduced states corresponding to = and fixed P^ can now be written 
as colour singlet Fock basis states on one link. A general state expanded in terms the orthonormal 
basis is then (summation on repeated indices implied) 

p=|n|,|n|+2,... 
p>0 
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•Tr{al(fci)at(fc2)...at(fcp)} |0) (31) 

where we set the winding number n = A + p+ -- - + f7 equal to zero. When the wavefunctions / are 
chosen appropriately, the states \'^) will describe an infinite discrete set of stable physical glueballs 
satisfying the mass shell condition. 

Using the mode expansion of the operators in Fock space, one may explicitly verify that the 
isomorphism under reduction holds. Some details of this calculation and its validity are given in 
Appendix The reduction rests on expanding about the M — classical solution. We have already 
argued that, at sufficiently small lattice spacing, the classical minimum is no longer at M — 0, and our 
quantisation is no longer valid. Note that this is not a statement about equivalence under reduction 
breaking down, but about our entire method, whether applied to the reduced or unreduced theory. 
Wilson's lattice gauge theory in the Lagrangian formulation is valid for any correlation length, but 
Eguchi-Kawai reduction breaks down at small lattice spacing [ p^ and various quenching and twisting 
prescriptions have been given in order to approach the critical surface in coupling constant space 
with a reduced theory. The dielectric formulation with expansion about M = does not need 
quenching/twisting, but we pay for this by having to work with an effective theory at rather large 
lattice spacing on the scaling trajectory. 

3.5 Symmetries 

The theory possesses several discrete symmetries. Charge conjugation induces the symmetry 

C : a^^i y- ^ al^ ji . (32) 

There are two orthogonal reflection symmetries Vi and V2 either of which may be used as 'parity'. 
If Vi : ^ —x^, we have 

Vi : a^i.,, - all,,, • (33) 

7^2- — > — x^, is more subtle in light-cone quantisation since it does not preserve the quantisation 
surface. Its explicit operation is known for sure only for a set of p free particles |2^ , which we call 

\ m' = l / 

The latter P2/ is nevertheless useful since it is often an approximate quantum number and its 
expectation value can be used to estimate V2 Given V2 and Vi we can determine whether 

spin is even or odd using the relation (—1)"^ — 'PiV2- If rotational symmetry has been restored 
in the theory, states of spin J ^ should form degenerate Vi doublets \ + J) ± \ — J) Q . We use 
"spectroscopic notation" | J'j^^'' to classify states. 
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3.6 The effective potential 

The basis written in is already diagonal in P+, and P+ commutes with . It remains 
to find the coefficient functions /, cyclically symmetric in their arguments, which diagonalise the 
Hamiltonian and hence the Lorentz invariant (mass)^ operator 2P+P^ with eigenvalue M^. 
Before proceeding we must choose the explicit form of the transverse effective potential. In this 
work, we will include in V[M] all Wilson loops and products of Wilson loops up to fourth order in 
link fields M . After Eguchi-Kawai reduction this limits one to the following form 

V[M] = /i^TrjMAfl'} + ^TrjMAftMAft} 

+^ Tr [MM'^MHi] + ^ Tr [mHi] Tr {mHi} (35) 

Note that the last term above, which might appear suppressed at large TV, is in fact non-zero only 
for 2-link Fock states (p = 2 in Eqn. (^l])). It is a term which 'pinches' the worldsheet swept out in 
light-cone time by the glueball flux ring given by the superposition of Wilson loops (|3^) . Near 
a critical surface in coupling constant space it is naively irrelevant, but at finite lattice spacing it 
turns out to be quite significant. In principle one should repeat calculations with higher orders of 
the field allowed in the action, to determine the magnitude of any changes. We will not be able to 
do this with our modest computing resources — there are many more allowed terms at 6th order — 
but will look at other quantities, such as expectation values of higher order Wilson loops, in order 
to assess effects of this truncation. (One can also of course compare with results obtained using only 
the quadratic term in (p5|)). Our goal is to find out if there is a trajectory in the parameter space of 
( p5| ) along which physical quantities scale and agree with the ELMC simulations ||] with reasonable 

The renormalisation of the quantum theory at fixed lattice spacing follows that of a 2D gauge 
theory with adjoint matter pO| |, involving only logarithmically divergent self-energy subtractions on 
the propagator and a principle value prescription for the light-cone Coulomb singularity. Explicitly 
one finds at large N 



fOO 

p+ ^ I dk 



fcTr{4(fc)aA(fc)} (36) 



P- = ^ I ^Tr{at(fc)aA(fc)} 



dk, 

2 Jo ^ 

[°° dki dk2 dks dkj 

+P / ,, , , , [5{ki + k2~k3~ ki)H2^2 

Jo 87ravfcifc2K3«:4 

+ S{ki +k2 + k3- ki) (Hi^3 + H3^l)] (37) 
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where 




1^3 • 



(38) 



(39) 



(40) 



In the N — oo hmit we do not distinguish U{N) from SU{N) and the gauge couphng g'^N is finite. 
The hght-cone Schrodinger equation defined by taking the matrix elements (Fock|2P+P^ j^E') = 
A/^ (Fockj^') allows one to write down an infinite set of successively coupled integral equations for 
the corresponding Fock components /. We do not bother to list them here since they are rather 
complicated, but refer the reader to Ref. pO[ | for equivalent expressions. To diagonalise we will 
employ a combination of variational ansatze for the /'s and numerical solution through discretisation 
the momenta k and truncation of the number of particles p. 

4 The glueball spectrum 
4.1 Analytic wavefunctions 

The behavior of the wavefunctions / in Eqn. (|3l|) when any one of the arguments vanishes 0] is 



and < /? < 1/2. Specifying also the number of nodes of / as a function of momenta, one can make 
a sensible ansatz. For Ai and A2 small there is very little mixing between Fock states of different 
number of link modes p | pO| . In this case a mass eigenstate | ^P) has predominantly a fixed p, the 
mass increasing with p. For a given p, the energy also tends to increase with the number of nodes 
in the wavefunction / due to the J (9-) J term in (p3), which is in fact a positive contribution. 



lim f\p...a{ki,k2, . . . ,kp) oc 



k" where 2/?tan(7r/3) 




(41) 



g^N 
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Thus one naively expects the lowest two eigenstates to be approximately 



dk 







/+i,_i(fc,P+-fc) Tr [al,{k)aU{P+ - k)} \0) (42) 



with the lowest state having a symmetric wavefunction /+i._i(fc,P^ — k), corresponding to 0^^, 
and first excited state having antisymmetric with one node, corresponding to . The 

next highest states are either a 4-link state with positive symmetric wavefunctions and 
or a symmetric 2-link state with having two nodes. In the glueball spectrum we 

identify the latter states as 0;^^ and 2++, respectively, although actual eigenstates are a mixture 
of these. The above naive considerations correctly predict the ordering of light glueballs in the 
spectrum, but differ in the details of the wavefunctions we find numerically later. 

In light of the cndpoint analysis (^), it is useful to consider wavefunctions / in Eqn. (|3^) with 
the simplified longitudinal form fcf and arbitrary transverse structure: a "transverse only" model. 
Thus we construct a basis of unit norm states of the form: 



■S {p+ - fc,) (fci • • • kpf (43) 



where ai G {—1, +1} and we set the winding number n ^ ct^ to zero. One can characterise these 
states with the "folding parameter" 7: 

P7 = The number of i such that ai = Oi+i. (44) 

Thus, 

p{\ — 7) = The number of i such that (Ji 7^ cr,;+i. (45) 

(Here, one should think of a 1-dimensional polymer model.) Roughly speaking, the transverse size 
of a typical state is a(p/2)^. Using 7, we can write down an exact expression for the expectation 
value of the Lorentz invariant (mass)^ operator: 

((Ti,...,(Tp|P^F''|cri,...,crp) 

- p(p(l + 2/3)-l). 

2 

= ^(52^ + 2(l-7)Ai+7A2) 
4a 

+ ^ {g^N (5 - 4 In 2) - (2Ai(l - 7) + A27) (4 In 2 - 1)) 

+0 (p) + O (/32) . (47) 
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For the two particle state |+1, —1) the expectation value of the invariant (mass)^ operator is, 

= ^ (.g'A^ + 2Ai + A2 + A3) 

+ — ig^N (3 - 2 In 2) - (2Ai + A2 + A,) (2 In 2 - 1)) 
a 

+0(/32) . (48) 

If one demands a positive-definite spectrum, these expressions place variational bounds on the al- 
lowed region in parameter space. First we must have positive /i^ else our quantisation procedure is no 
longer valid and the spectrum becomes unbounded from below. For the p = 2 case the expectation 
value of the invariant (mass)^ becomes tachyonic when 

2Ai + A2 + A3 < -g^N (1 + 8/3) + O (/J^) , (49) 

while from (|4^) we see that many particle states become tachyonic at 

2(l-7)Ai+7A2 < -ff27V(l-f8/3)-HO(/3') . (50) 

4.2 Numerical results 

In our numerical solutions we restrict the number of link fields p in our basis states (|3l| ) and discretise 
momenta by demanding antiperiodicity of the fields in ^ + L. This is a high light-cone 
energy cut-off since, unlike in the transverse direction, light-cone momentum goes like the inverse 
of energy, (we use the method of to handle the IR divergence). For fixed integer valued cut-off 
K = LP+/(27r) momenta are labeled by odd half integers Km ~ Kk„i/P'^, K — J2m '^"i- C)n the 
computer, we generate a basis of states and calculate matrix elements of the Hamiltonian using the 
hybrid wavefunction/DLCQ matrix elements discussed in Appendix]^. The nonzero matrix elements 
are then stored in a file (the matrix is, in fact, quite sparse). To diagonalise the Hamiltonian P^, 
we employ a standard Lanczos algorithm. 

In order to fix the coupling constants in the effective potential, we perform a least fit to 
Teper's ELMC large N extrapolated spectrum]^ As we shall later show, the mass in units of the 
coupling = ii^a/{g'^N) is a measure of the lattice spacing while the other terms in our effective 
potential Ai, A2, and A3 must be found from the fitting procedure. We will also determine g'^N/a 
based on a fit, which we check self-consistently with measurements of the string tension. 

In order to minimise errors associated with our truncation in p and K, we take the spectra for 

{K,p truncation) G {(10, 4), (10, 6), (10, 8), (11, 6), (12, 6), (13, 6), (14, 6)} (51) 

^For the ELMC data, we assumed rotational symmetry and averaged over parity doubles. 
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Figure 1: Difference between tfie unextrapolated and extrapolated eigenvalues A = 

=i4,p<6 -^"^extrapolated ^'^^ Spectrum in Fig. ^. Also shown is the function 0.05M^ 

which is the assumed size for our error bars. 
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Figure 2: Scaling trajectory in parameter space. We plot (a) Xi/{g'^N) and (b) \2/{g'^N) vs mass 
m^. The solid line is the data in Table |l]; the dashed line is from a second trajectory discussed in 
Section 4.3. 



and extrapolate to the continuum, K,p oo, using a least squares fit to the function 



C2 e 



cp 



(52) 



where c is given by our estimate in Appendix |^ and the 1 /K term is consistent with the leading 
finite K error as discussed in Appendix Since we cannot measure \J^\ directly, we only classified 
states according to Vi and C during the fitting process. In practice we found that the fitting function 
( |52|) worked quite well for all but the highest 3 or 4 states where the exponential convergence in 
p-truncation was not yet evident.^ We estimate our one-sigma errors from finite K and p-truncation 
to be roughly O.OSM^; see Fig. ^ 

As we fit the various couplings as a function of m^, we find a narrow strip in parameter space 

^However, looking at the K = 10, p < 8 sector, even these highest states had less than 15% 8-particle content. 
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^1 

^2 AT 


^2 

^2 AT 


y 

ac 


x' 


c 


0. 


-0.00565 


-0.305 


4.38 


30. 


-1.43 


0.00394 


-0.035 


-0.366 


4.11 


25.9 


-1.29 


0.0158 


-0.0677 


-0.433 


3.84 


25.3 


-1.2 


0.036 


-0.123 


-0.496 


3.65 


23.1 


-1.1 


0.065 


-U.zOz 


-0.55 


0.4/ 


23.5 


-i.Oi 


0.104 


-0.297 


-0.602 


3.28 


23.6 


-0.94 


0.153 


-0.427 


-0.626 


3.1 


24.9 


-0.884 


0.214 


-0.58 


-0.658 


2.97 


25.7 


-0.832 


0.291 


-0.778 


-0.657 


2.84 


28.6 


-0.792 


0.384 


-1.05 


-0.508 


2.67 


35.2 


-0.802 


0.5 













Table 1: Scaling trajectory for the extrapolated spectrum. Here we tabulate couplings vs for a 
best fit to the ELMC spectrum. We find that \^l{g'^N) is always quite large (> 100). We could 
not find a satisfactory minimum for the rn? = 0.5 datum. 

where we obtain good agreement with the ELMC spectrum; see Table |l| and Fig. |. As we shall 
see, moving along this strip corresponds to changing the lattice spacing a. The strip, where has 
a local minimum, disappears when m? is sufficiently large, indicating that for large enough lattice 
spacing our truncation of the effective potential is no longer a good approximation. 

In Fig. ^, we have plotted a typical spectrum along the scaling trajectory together with the 
ELMC results. For graphing purposes, we assigned \ J\ to our spectrum based on a best fit to 
Teper's results. As a partial check of this assignment, we calculated the expectation value of 7^2/, 
Eqn. (Q), for our spectrum {K = 10, p < 6 truncation). The subsequent determination of (— 1)"^ 
agreed with the spectrum fit for most of the low-lying glueballs. However, there was disagreement in 
the assignment of (—1)'^ between the lowest ^ , first excited 2 and lowest 1 ^ states in Fig. ||. 
However, since these three levels are almost degenerate, this discrepancy is not important. Of more 
importance is the \ J'\^ sector where our measurements were ambiguous (with (5'|7^2/ \^) ~ 0) 
except for the lowest 2^ state where we obtained the opposite of the desired result for (— 1)"^. 

Although the overall fit with the conventional lattice results is quite good, we see two deficiencies 
of our spectrum that cannot be attributed to K or p-truncation errors. First, we see that the lowest 
0~~ state is too low in energy. Second, we see that the lowest parity doublet 2*+ is not quite 
degenerate. We believe that these discrepancies must be due to our truncation of the effective 
potential. Finally, we have made no prediction for the lowest 1^^ state since it lies too high in 
the spectrum; there may be a additional 0+^ state of lower energy which prevents us from 

labeling the 1++ unambiguously. 
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Figure 3: A fit of our LCTL extrapolated results against Teper's ELMC large N extrapolated 
spectrum 1^. The eigenvalues are shown in units of (Teper's) string tension for various 
The parameters are from the — 0.065 row of Table ^. One finds similar spectra along the entire 
scaling trajectory in coupling constant space. 




0.5 



Figure 4: Mass — iJ?a/{g'^N) such that the lowest eigenvalue is zero vs \i/{g'^N) and 
X2/{g^N) for p < 6, K = 14, = IQQg^N. Below this surface the spectrum is tachyonic and below 
771^ = our quantisation breaks down. The dark line is the scaling trajectory. 
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M2 

a 

36.2 

67.3 

27.3 

50.5 

73.5 

109. 

86.3 

38.8 

82.4 

51.8 

94.9 

67.1 

102. 

67. 

124. 

92.2 

79.4 

84.4 



0.036 
0.046 
0.051 
0.909 
0.06 
0.02 




0.86 
0.022 





0.859 
0.054 







0.047 




0.472 
0.359 
0.405 
0.075 
0.681 
0.16 
0.142 
0.244 
0.09 
0.113 
0.885 
0.16 
0.082 
0.609 
0.909 
0.701 
0.549 
0.398 
0.895 



0.372 

0.528 

0.332 

0.011 

0.153 

0.741 

0.026 

0.671 

0.033 

0.096 

0.014 

0.036 

0.05 

0.132 

0.003 

0.126 

0.156 

0.447 

0.007 



0.019 

0.019 

0.026 

0.002 

0.037 

0.014 

0.026 

0.01 

0.006 

0.078 

0.027 

0.011 

0.003 

0.061 

0.037 

0.027 

0.022 

0.013 

0.034 



0.084 
0.024 
0.145 
0.002 
0.065 
0.048 
0.801 
0.059 
0.011 
0.512 
0.074 
0.793 
0.005 
0.136 
0.051 
0.142 
0.269 
0.081 
0.063 



0.017 
0.023 
0.04 


0.003 
0.018 
0.005 
0.016 


0.179 


0.001 
0.001 
0.007 


0.004 
0.003 
0.014 




Table 2: The transverse structure of the spectrum for K = 14, p < 6. The first column is our 
assignment of \J^\^^'^, the second column is the extrapolated (mass)^ from Fig.|3[ and the remaining 
columns show the probability for a state to have a given transverse structure. The coupling constants 
are from the = 0.065 row of Table § 

The numerical bound for absence of tachyons is shown in Fig. ^ as a zero-mass surface. The 
surface agrees rather well with the analytic bound provided by the "transverse only" model ( |50| ) 
with 7 taken over the full range < 7 < 1. As the transverse lattice spacing vanishes the mass 
gap should vanish in lattice units. The fixed point for this, which we believe lies somewhere at 
negative m^, should lie on the zero-mass surface, but is inaccessible to us in the dielectric regime 

> 0. Nevertheless the scaling trajectory should gradually approach the zero-mass surface if it is 
to eventually encounter the fixed point. We see this happening already in Fig. ^. 

Looking at the transverse structure in Table ^, one can see several important properties of the 
spectrum. First we notice that the lowest O"*""*" state has very little two particle content. This is the 
main effect of the (Tr {M'''M})^ interaction when the associated coupling constant A3 is large. It 
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essentially removes the two particle sector from the wavefunction of this state. Second, we notice that 
there is a relatively large mixing between the various sectors, including mixing of sectors of different 
particle number at the 10-20% level. This is large enough to restore some rotational invariance but 
not so large that it contradicts our truncation of the effective potential. If we look at parity doublets, 
we see that each wavefunction is quite dissimilar from its partner, including very different particle 
content. This is what one would expect since full parity is dynamically complicated in light-cone 
coordinates. Finally, one would expect from Regge trajectories that the lowest 2++ state would have 
a large two particle content while the first excited state should not. This helps us label the 
correct \ J\ for these two nearly degenerate states. In fact, large mixing between these states may 
explain why we do so poorly with the lowest 2^'^'^ parity doublet. 

4.3 Is our spectrum accidental? 

Since the scaling trajectory in Table ^ is a local minimum of in coupling constant space, one 
might worry that our agreement with the ELMC data is accidental. Since we are not always able to 
assign \ J'\ to our states, it is possible that another local minimum in parameter space represents the 
'correct' scaling trajectory. With this in mind, we have conducted a thorough search of parameter 
space and have found two other trajectories. 

The first occurs for A2 ~ 0.5g^N and Ai < 0. However, the of the fits are significantly 
worse. Moreover, it goes against our theoretical prejudice that the scaling trajectory should be in 
the vicinity of, and track towards, the zero-mass surface (fig.(^)). 

The second occurs quite close to our scaling trajectory. We have included Ai and A2 for this 
trajectory as dashed lines in Fig. |2[ Indeed, the of the fits are as good or better than that of 
our chosen trajectory. We find that the principle difference between the two trajectories is that the 
2++ and 0'*'"'' levels are switched — recall that we are unable to distinguish these levels by exact 
symmetries alone. Consistent with the Regge trajectory picture, we believe that two particle states 
should form the lowest Regge trajectory and hence the lowest 2++ state should contain mainly two 
particles. Thus, we have chosen the scaling trajectory where the lowest 2++ state has mainly two 
particle content. Additionally, we find that the ratio (g^cr)LCTL/(3^c)ELMC in Fig. ^ is not as good 
for this other trajectory (dashed line), nor does it appear to track towards the zero-mass surface. 

5 The String Tension. 

To measure the string tension in the direction (before Eguchi-Kawai reduction) consider a lattice 
with n transverse links and periodic boundary conditions. Constructing a basis of Polyakov loops, 
"winding modes," that wind once around this lattice, one may extract from the lowest eigenvalue 
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Af^ vs n the lattice string tension aa — AM„/An. Because of Eguchi-Kawai reduction, this is 
equivalent to using Polyakov loops of winding number n on the single- link periodic lattice (Eqn. ( |3l| ) 
with n 7^ 0). 

5.1 Analytic results 

The behavior of Eqn. (^l|) is equally true of winding modes n ^ 0. Thus, it is useful to extend 



the "transverse only" model of Section 4.1 to the case of nonzero winding number n. In particular, 
consider thep — n particle state |-|-1, +1, • • • + 1) and thep = n+2 particle state | — 1, +1, +1, • • • + 1). 
The correct normalisation of these states is 

1 = {+1, + + ■■■ + !) =pCl^.+i....,+i (53) 

1 = (-1,+1,+1,--- + 1| -1,+1,+1,--- + 1) 1^+1^+1^ +1 . (54) 

Using these two states as a basis, we obtain a 2 x 2 mass squared matrix (n > 1) 

^ / 7rn(n - 1) (g^AT + A2) 4nV^rTT(Ai + A2) \^n(R\ (^^\ 

ia\ 4nV^rTT(Ai + A2) Tr{n + l){n + 2) {g^N + X^) + Awin + l)Xi J ^ ' ^^^> 

The diagonal elements come from Eqn. (p6|), setting 7=1 and 7 = {n + l)/{n + 2), respectively. 

The off-diagonal matrix element is, more precisely, 

(-f 1, +1, • • • + llPuP^I - 1, +1, +1, • • • + 1) 



(Ai + A2) r(4/3 -I- 1) r(/3 + 1) Vnr ((2/3 + l)n) T ((2/3 + l)(n + 2)) 
2™ r(3/3 + I) r (2/3 1)2 r (2/3(n -f 1) + n) 
An example of the lowest eigenvalue of Eqn. ( |55| ) vs n is plotted in Fig. |[ For large n one obtains 
a linear spectrum M (A2 + g^N") /a and the vanishing lattice string tension condition A2 = 

-g^N + 0{(3). 

5.2 Numerical results 

Numerical calculation of the winding modes proceeds in the same manner as our calculation of 
the spectrum except that we use a basis of n 7^ states. String theory arguments indicate that 
oscillations of the winding mode transverse to itself yield a form Ad'^ — o^^naY ~ (ttt/S, the constant 
correction being due to Casimir energy p3|] . Fig. ^ shows a typical vs n plot for winding modes. 
Consistent with the expectations from string theory, we have fit the winding mode spectrum to a 
curve of the form 

^ = a-on- + , (57) 
a a 

which we expect to be appropriate for QCD at sufficiently large n. As illustrated by Fig. |^, there is 
a good fit to the quadratic. In the same manner, we determine ( [5^ ) all along the scaling trajectory 
with the result is shown in Table ||. 
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Figure 5: Lowest eigenvalue vs n for winding modes. Here, K ~ 10.5 or 11, p < ?i + 4, and the 
couplings are taken from the m} — 0.065 row of Tablejl]. Also shown is a numerical fit to 1.92n^ — 3.5. 
The dashed line is from the variational estimate of Eqn. (p5|). There is a subtlety associated with 
the n = 1 point, which is described in the text. 
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0. 


1.48 


-2.77 


0.00394 


1.56 


-2.95 


0.0158 


1.65 


-3.08 


0.036 


1.77 


-3.29 


0.065 


1.92 


-3.5 


0.104 


2.09 


-3.68 


0.153 


2.3 


-3.82 


0.214 


2.55 


-4.06 


0.291 


2.86 


-4.28 


0.384 


3.29 


-4.21 



Table 3: String tension measurements along the scaling trajectory. The coupling constants are 
taken from Table |^. Here, the lowest winding mode eigenvalue M^/cr has been fit to a^av? + b/a 
for n = 2,3,4,5. 
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Figure 6: The ratio (g^cr)LCTL/(5^cr)ELMC- The sohd hne is for our scaling trajectory while the 



dashed line is from the second trajectory discussed in Section 4.3 



The constant term b/a represents the leading finite transverse size correction to the string tension, 
and is supposed to be universal. The fact that it varies strongly along our 'scaling trajectory', shows 
that for this quantity there are still large finite a errors. Also we have not extrapolated in K or p 
truncation. There is a subtlety associated with the n = 1, one particle matrix element 
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Tr{ai(X)}P^P^Tr{al(if)} |0) 



(58) 



The most consistent definition is to assume that this matrix element is infinite due to the diagonal 
part of the J(9_) ^ J interaction. In this case, M^^^ lies somewhat above the extrapolation of 
the solid line in Fig. ^ down to n — 1, since the eigenfunction contains only 3, 5, 7, .. .-particle 
contributions. (If one instead defines ( |58| ) to be equal to /i^, M^^i falls right along the extrapolation 
of the quadratic fit.) We have not included the n = I datum in our determination of the string 
tension. 

If we are careful to distinguish between our measurement of g^N and a versus g'^N and a used 
in the ELMC calculations, we can form the ratio 

(g''^)LCTL ^(^] (g'^)LCTL / "''^gcTL (59) 

The first term is 0.1974(12) according to Teper's large N extrapolated measurement of ^/a/{g'^N) 
and the remaining two terms are from Tables ^ and ^, respectively. This ratio measures our error in 
measuring the string tension plus any discrepancy between g'^N as used by us and that used in the 
ELMC calculation. From Fig. |[ we see that the two couplings differ by 5%. This implies a rather 
large discrepency since we have used a non-standard definition of g^ . 

In general, we find that the surface in our parameter space on which M"^ vanishes (in units of 
g'^N/a) , Fig. H, always lies above the surface of vanishing lattice string tension a^cr; that is, the 
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Figure 7: The lattice spacing in units of the string tension along the scaling trajectory, Table | 
Here we plot a^a vs rn^. 



spectrum always becomes tachyonic before the lattice string tension vanishes. Any "critical point" 
must occur in the intersection of the surface of vanishing lattice string tension and vanishing mass 
gap. In fact, we find no such point anywhere in the > coupling constant space, in agreement 
with the colour dielectric arguments of Section ^ Instead we have a scaling trajectory at finite 
lattice spacing, as we now show. 

If we assume that our cr is equal to the ELMC value of cr, we can determine the lattice spacing a in 
units of the string tension along the scaling trajectory. As mentioned previously, this demonstrates 
that the mass — fi^a/{g'^N) determines the lattice spacing and that the continuum limit occurs 
at /Lt^ < 0; see Fig. 0. 

6 Density of States and T^. 

Another quantity of interest is the deconfinement temperature T^, which on general grounds one 
expects to exist in a pure gauge theory in greater than 1 + 1 dimensions . In the Euclidean finite 
temperature formalism, this is signalled by a winding mode in the Euclidean time direction which 
becomes tachyonic. An estimate of the inverse circumference T = Tc at which this occurs comes 
from the formula [p5| for the lowest mass eigenstate of a long closed harmonic string of winding 
number one: 



2 - 2) 

M 

where D is the number of spacetime dimensions. It becomes tachyon at 

/ 3a 



t:{D-2) 



(60) 



(61) 
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The result Tc w 0.98y/a for D — 3 is close to, but different from, that predicted by the large-iV 
extrapolation of the ELMC results |||:[| = 0.89(3) V^. 

Tc may also be obtained from the Hagedorn behaviour of the asymptotic density of mass 
eigenstates 

p{M) - eyi^{M/Tc) (62) 

The canonical partition function diverges for T > Tc. If a > (D + l)/2 it is a phase transition, 
beyond which the canonical and microcanonical ensembles are inequivalent. If a < (_D + l)/2 the 
ensembles are equivalent and Tc represents a limiting temperature — the free energy diverges at Tc- 
The light-cone quantisation of a free relativistic string with harmonic oscillations in the transverse 
direction (Nambu-Goto Bosonic string) was first performed by Goldstone et alii [|2^, who found a 
linearly rising discrete spectrum of squared masses 

- 27rcrn , n£Z (63) 

very similar to the hadronic Regge trajectories observed in nature. Each mass level has a degeneracy 
increasing with n, such that the density of states p per unit mass yields Eqn. ( |6l| ) once more. Thus 
large- QCD seems to have slightly larger level density than this Bosonic string. 

In the next section, we will use the density p{M) in the LCTL formulation to understand an- 
alytically this difference between Tc in the large N gauge theory and the Bosonic string. We will 



discuss numerical estimates of Tc from the Hagedorn behavior in Section 3.2. 
6.1 Analytic Regge Trajectories. 

High mass glueballs are not quite described by the Bosonic string model. One piece of evidence 
for this is the discrepancy in Tc noted above. By construction, the only degrees of freedom, in 
addition to the centre of mass, allowed in the Bosonic string model leading to ( |6l| ) are transverse 
oscillations. A 'QCD string' might also be expected to have longitudinal degrees of freedom. In 
fact it is well-known that allowing a relativistic string to oscillate only in transverse directions leads 
to consistency problems, for example with Lorentz invariance, unless D = 26. A correct treatment 
when 13 7^ 26 must allow longitudinal oscillations in an essential way Recent attempts, by one 
of the authors, to study directly the light-cone quantisation of such 'non-critical' string models have 
led to largely inconclusive results however ||3^ . The Wilson loops of our light-cone lattice dielectric 
theory clearly exhibit oscillations in the both the longitudinal and transverse directions. We 
propose that the former accounts for the slight shift in Tc- 

^The numerical ELMC result is in turn in agreement with the analytic large- Af result of the Torino group for Tc 
based also on Euclidean lattice methods 
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Figure 8: Coulomb exchange diagram and its self-energy contraction. Thick lines are propagating 
particles, wavey line the non-propagating mode. 



It is instructive to derive the Hagedorn behaviour in the LCTL approach and show how the extra 
degrees of freedom necessary to lower from its Bosonic string value might arise. The exponential 
density of asymptotically high mass eigenfunctions of ( p^ derives from wavefunctions containing 
typically many links p, i.e. long strings. It is useful to consider two extremes of this large p problem, 
according to the complexity in the longitudinal dependence of the wavefunctions / on Bjorken 
X — k/P^ , conveniently measured by the typical number of nodes say. We split P^ = P^ + Pf 
into transverse and longitudinal pieces, where from Eqn. (e3) 



There are two problems which face us now. Firstly, we expect that more terms must be added to 
y[M] as we examine levels of higher mass in order to restore Lorentz invariance for them. Secondly, 
the wavefunction will not in general split into a product of longitudinal and transverse pieces. There 
are two extremes when one might expect a degree of universality to help. 

Consider first the extreme when p takes its maximum value, for example when each link variable 
assumes the minimum x = 1/K allowed by the DLCQ cut-off. The solutions of P^ in this limit 
have already been studied by Klebanov and Susskind |3|] (see also Ref. in the approximation 
that processes which lower the number of particles are neglected.^ They found that for a wide range 
of parameters in the potential V the spectrum was precisely (^3|) with the correct degeneracy to 
produce (^l|). Pf^ vanishes in the approximation, so one essentially diagonalises P^^ . They also 
argued that allowing some particle annihilation and creation would merely renormalise cr to a first 
approximation. This would mean that, as far as the asymptotically high spectrum is concerned, 
the purely transverse degrees of freedom of glueballs in the LCTL formulation are just those of the 
Bosonic string. 

In the previous extreme the longitudinal dynamics were frozen by construction. Let us now 

*We will not regurgitate the somewhat technical analysis, but refer the interested reader to the original literature. 




(65) 



(64) 
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unfreeze them by allowing oscillations in the longitudinal direction, with non-trivial dependence of 
the wavefunctions on Bjorken x, and go to the other extreme. The high spectrum for Pj^ typically 
contains wavefunctions f\p...aixi, . . . , Xp) which oscillate rapidly as a function of each Xi. In this case 
we can make the asymptotic approximation suggested by 't Hooft for fundamental representation 
fermions in two-dimensional light-cone large-iV QCD; in fact our initial analysis is only a slight 
modification of the corresponding one for adjoint representation fermions due to Kutasov 
For such rapidly oscillating wavefunctions all of the terms in tend to average to zero. The same is 
true of the terms in Pj~ , except for the Coulomb exchange term which is singular at zero momentum 
transfer (fc2 ^3 in Eqn. (|3^)), and the creation of small-x particles from the same instantaneous 
mode (Eqn. ( p9| ) at fci -|- ^2 = or fc2 + fcs = 0). Keeping only the Coulomb exchange term 
for the moment, the Pf eigenvalue problem is dominated by the following part of the light-cone 
Schrodinger equation for the wavcfunction components of a p-link state 

M JXp...a[Xl, ■ ■ ■ ,Xp) = / -^{jXp...cr[Xl, ■ ■ ■ ,Xp) 

™ J-oo y 

— fxp...a{xi + y,X2 — y, ■ ■ ■ , Xp)} + cyclic permutations (66) 

This is obtained by projecting the Coulomb part of Pj^ onto individual Fock states. The Coulomb 
kernel has been simplified and the limits of the integral extended to the whole real line under the 
assumption that the dominant contribution comes from y w 0. The right hand side of ( |66| ) is 
represented by the diagrams of light-cone perturbation theory in Fig. ^. We must decide what 
to do about the transverse shape of the Wilson loop, represented by the indices (A, p, . . . , cr). In 
principle there will be a fine structure spectrum of transverse shapes on top of the large mass 
eigenvalue provided by Eqn. (|66|). We will make the approximation that this can be represented by 
an overall transverse degeneracy factor, which for simplicity we will take as the number of transverse 
configurations of a p-link loop ^ 2^. This is an overcounting, obviously including lattice artifacts, 
and a more sophisticated treatment is no doubt possible. With the random distribution of transverse 
configurations, we can suppress the indices (A, p, . . . , cr) for the following; however, the number of 
links p is remains significant. Solutions to (^6|) must satisfy the additional boundary conditions 

/(xi, a;2, . . . , Xp) = if any = , me {!,..., p}. (67) 

If we look for solutions symmetric under a cyclic permutation of the x^'s, a sufficient condition 
for (^) is that /(O, X2, . . . , Xp) = 0. The solutions are essentially those given by Kutasov |3^] with 
minor modifications for Boson statistics. We just display the first couple: 

f{xi,X2) — sin7rrtia;i ; — 2g'^NTrni/a ; ni odd integer (68) 
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Figure 9: Self-energy corrections to the non-propagating mode. These diagrams represent dom- 
inant particle non-conserving processes when small momentum fraction y is flowing through the 
initial line. 

f{xi,X2,X3,X4) = sin(7rni(a;i + 2:2)) sin (7m2 (x2 +2:3)) - sin (TrniXi) sin (7^2X3) 

-|- sin (7:711X2) sin (7771,22^4) + ^ "•2) ; (69) 
— 2g'^NTT{Tii +n2)/a ; rii > n2 ; Ui odd integers. 



The general spectrum is 



= ^^-J^{ni +n2-\ h np/2) ; Hi odd. (70) 



In general the high mass states are a mixture of different numbers of particles p, so Eqn. ( |66| ) 
cannot be the whole story. We must also take into account pair creation effects. The most important 
effect is likely to be the A+ self-energy processes Fig. (|), which are singular when small momentum 
fraction y is flowing through the instantaneous line. The net effect of such "vacuum polarisation" 
is to give a renormalisation of the Coulomb coupling constant g^Njira and hence the mass scale 
setting for Eqn. ([70|). In fact such a renormalisation must occur since the bare expression ( [70| ) is 
obviously a dependent; the Coulomb exchange process cannot produce any motion on the transverse 
lattice. Through pair creation and annihilation of links, transverse motion becomes possible. We 
will not need to determine the renormalisation for the present application, since the spectrum is still 
of the string theory form ( |63| ) and our ignorance of the functional dependence of the mass scale on 
the couplings is parameterised by a. 



From the form (70) we may calculate the density of states using standard methods of statistical 
mechanics. The generating function 

G(zi;) =^d„zi;" (71) 

n 

defines the number of states dn of Af ^ — 27zan. For large n we need to find the behaviour w ^ 1. 
Including the transverse degeneracy factor 2^ we have explicitly 

oo 

logG(u;) = log [| (1 + 4^") 
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where Li2 is the dilogarithm function. Then 



(72) 



dn — - — / dw -"T (73) 



1 /..,G(«^) 



where C encircles the origin. The integral ( [73| ) has a saddle at logw = ^J^^LA2{—^)/n giving 

(74) 



exp 



2V-nLi2(-4) 



exp 



The Hagedorn temperature Tc — O.SI^/ct should be regarded as a lower bound given the probable 
overestimate of the degeneracy of transverse configurations. The important point to note is that the 
density of states would rise exponentially as a result of longitudinal excitations alone, even if we did 
not include any degeneracy factor for transverse structure. This is strong evidence that Tc can be 
lower in large- A'^ gauge theory than in free Bosonic string theory as a result of non-trivial longitudinal 
degrees of freedom. Our estimate, therefore, is that it lies in the range 0.81^/a < Tc < 0.98-^^, 
consistent with the ELMC result. 

6.2 Numerical estimates 

Now we turn our attention to a numerical estimation of the density of states. For this calculation, 
we have only specified the box size K to truncate the basis. It is essential to demonstrate that the 
spectrum is sufficiently converged in K, thus we only fit to the states between 0.5 < log(i) < 5.5 for 
the K — 7 data; see Fig. A numerical fit gives, 

log(t) =3.99+ ^-5.86 log (il). (75) 

Since the density of states is p{M) — dt/dM, we find Tc — 0.78y/a with an estimated erroij^ of 
at least 10%. Due to the large error, this result is both compatible with the Euclidean lattice 
result and the analytic bound of the previous subsection which assumed random transverse 

configurations. The fact that the power correction a ~ 4.86 in the density (^) is much larger than 
{D + 1)/2 ~ 2 means that the thermodynamic free energy remains finite at Tc, marking a true phase 
transition (presumably deconfinement) rather than a limiting temperature. 

^This result could be easily improved upon since the K = 7 truncation contains only 4117 states. 
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Figure 10: Int vs the mass of the t-th eigenvalue Mt in units of the string tension. Here we 
have appUed a cutoff in K only. Also shown is a least squares fit to the K = 7 data in the range 
0.5 < log(i) < 5.5. We have chosen the couplings from the = 0.065 line of Table |^. 
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7 Conclusions 

We performed light-cone quantisation of the colour dielectric lattice gauge theory on a transverse 
lattice, for 2 + 1 dimensions and large N. The glueball solutions of the theory were investigated as 
a function of couplings in an effective potential T^[M] on the transverse lattice, which contained all 
possible gauge-invariant terms up to 4th order in M; 4 couplings in total (plus an overall mass scale). 
We found an approximate scaling trajectory — a one-parameter set of couplings for which physical 
results are almost unchanged — along which our glueball spectrum agreed quite well with results 
from Euclidean lattice Monte Carlo simulations; some 18 or 19 energy levels. The surface of vanishing 
mass in lattice units was also found in the region of parameter space where our quantization is valid. 
Our scaling trajectory correctly moves towards this surface as we decrease the lattice spacing. It 
presumably differs from the exact scaling trajectory by higher order terms in y [Af] that we neglected 
and which represent at this stage a systematic error. However, these corrections are seen to become 
less important as one moves along the trajectory toward smaller lattice spacing. We also bounded and 
estimated the deconfinement temperature from the density of asymptotically high masses, suggesting 
that longitudinal string-like oscillations are responsible for deviations from the Nambu-Goto string 
model. Our methods in addition yield the light-cone wavefunctions for every glueball in terms of 
the dielectric link variables, a more detailed analysis of which will be presented in future work. We 
applied new numerical methods for the solution of light-cone Hamiltonians which should be useful 
for a range of other problems. Below we mention a couple of important questions arising from our 
analysis and discuss improvements and extensions to this work that we are carrying out. 

7.1 Large- phase transitions? 

In lattice gauge theories with a simple action, such as the Wilson action, there is a strong crossover 
behaviour between strong and weak coupling which probably turns into a phase transition in the 
large N limit. This large N phase transition is well established in exactly solvable models [Q and 
is a lattice artifact, id est it occurs even on a single link. (Remember there are an infinite number 
of degrees of freedom per link at large N.) We have been tacitly assuming throughout this work 
that the infrared stable RG trajectory does not encounter such a large- iV obstacle between the large 
lattice spacing regime (where the dielectric approach is valid) and the infra-red fixed point. Any 
truncation of the effective potential almost certainly will encounter a large- iV transition. This 

is simply illustrated by the following exactly solvable one-link theory |3q| 




(76) 
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where W{x) — g2X^ +53^;'^ + • • • is a generic polynomial potential, with all > (this condition is 
convenient but not essential). Decomposing M = DV into a real diagonal matrix D and unitary 
matrices U and V, the solution can be expressed entirely in terms of the distribution of eigenvalues 
p(A) dX, which is determined by a saddle point analysis at large N. Here, A > is an eigenvalue of 
M^M. For > the support of p(A) is [0, 1], under suitable normalisation of M, and p(A) ^ 1/VA 
as A 0. For sufRciently negative /i^ the support becomes [A, 1], A > 0, with p(A) ^ \/X — A as 
A — !■ j4. This change in the analytic structure causes a third order phase transition in Z at a certain 
(negative) value of fi^ , in direct analogy with the situation for Wilson link variables [Q . In order to 
avoid the phase transition, one must ensure the support exists down to A = 0. On the other hand, 
according to Section |^, the effective potential should force M 1 as /i^ ^ — oo. These two 
requirements can only be met by carefully tuning a non-polynomial effective potential. The infrared 
stable trajectory would be of this sort. The density p would begin to vanish everywhere, except 
near A = 1, only as the fixed point itself is reached. Any finite-order truncation of the effective 
potential would not leave a navigable route to the fixed point, although it may allow a trajectory at 
intermediate lattice spacing close to the infrared stable trajectory. It is this latter possibility that 
we have sought to establish in this paper. 

7.2 Zero modes? 

There are two sources of fc+ = modes in the transverse lattice theory, both of which we have not 
used. Firstly the light-cone gauge A_ = is not compatible with removal of the region fc+ = 0, 
although the similar gauge (9_A_ = is. The difference is an extra quantum mechanical degree of 
freedom which will lead to a non-trivial vacuum structure in the type of theories considered here 
[ p6| . Very little is known about the quantitative effects on the spectrum. Some of the effects may 
just be absorbed into our effective potential. This dynamical zero mode certainly deserve further 
attention. 

Secondly, there is the — part of M, which we can ignore in the dielectric regime. If 
we decompose M — HU into a 'radial' Hermitian matrix H and an 'angular' unitary matrix U, we 
expect that as the mass squared of the Af -quanta turns negative, H gets a VEV Hq. If iJ = Hq + H, 
as the lattice spacing vanishes Hq freezes while H becomes infinitely massive and decouples, leaving 
U which is nothing but the usual Wilson link variable. This is just another way of phrasing the 
arguments of Section ^ Evidently the Hamiltonian approach becomes very complicated once outside 
of the dielectric regime because of this zero mode and the neccesity of decomposition into radial and 
angular parts. The fact that we see a slight worsening of our results at very small /i^ may indicate 
that we are becoming sensitive to the onset of the above phenomenon. Howvever we cannot rule out 
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that these fluctuations are due to some subtlety in the K and/or p extrapolations near /i^ = 0. 
7.3 Improvements and extensions 

There are two points of view one could take with regard to fixing the couplings in the effective 
potential y[Af]. If one is mostly interested in determining a hadron light-cone wavefunction, an 
acceptable phenomenological procedure would be to tune y [M] until it reproduces the mass spectrum 
from experiment or Euclidean lattice methods. The hadronic structure functions et cetera would 
then be predictions. We will not discuss this point of view further. Another approach is to attempt 
to determine the scahng trajectory a priori without reference to other data by examining Lorcntz 
invariance. In this paper our only measure for this was the multiplet structure of the glueball 
spectrum. Although in principle one could fix V^[Af] by demanding the correct degeneracies, in 
practice we found that this was not really tenable since we are working at finite lattice spacing. 
Indeed we found it difficult to get the lowest parity doublet to be degenerate. In addition, it takes 
no account of rotational invariance in the spin states, which form most of light spectrum. Another 
method is to measure the heavy source potential in the spatial directions. In Ref. this method 
was studied for the same light-cone Hamiltonian problem, but with arbitrarily set to zero. 

Nevertheless, for distances greater than two lattice spacings the potential appeared remarkably 
rotationally invariant in {x^,x^). However, over a typical hadron size (of order two or three lattice 
spacings according to our analysis) the results were much less satisfactory. Also, while it can test 
the scalar part of the interaction, it is not sensitive to the full spin dependence of the theory. 

A better measure of the restoration of Lorentz invariance is provided by the relativistic dispersion 
relation of each glueball wavefunction which allows one to keep track of Lorentz invariance for all 
states, whether of zero or non-zero spin. To test whether the dispersion relation A'P = 2P+P^ — 
{P^Y holds for each glueball requires us to work at non-zero transverse momentum Pi. This can 
be done within the context of the Eguchi-Kawai reduced formulation by adding appropriate phase 
factors to the Hamiltonian matrix elements. It does require a certain amount of code rewriting and 
we plan to include this measurement in future work. It should enable us to fix an optimal effective 
potential for symmetry restoration and thus predict the approximate scaling trajectory without the 
need for a fit to an existing spectrum. 

Although we have not discussed them in any detail in this paper, perhaps the most interesting 
output of our calculations are the light-cone wavefunctions. These may be used to calculate not only 
distribution functions at a given low normalisation point for the transverse momentum cutoff, but 
also their non-perturbative evolution as one moves along the scaling trajectory. A detailed analysis 
will be given in a future publication once the trajectory has been better verified. 
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There is the possibiUty of performing some sort of analytic or semi-analytic renormalisation group 
procedure. As we have emphasised, we cannot use small lattice spacing as a parameter to do this 
in the dielectric regime. However, we do have dynamical information which suggests that mixing of 
Fock sectors with different numbers of particles can be treated as a perturbation, and we can also use 
perturbation in small if necessary. In this paper we have already performed some rudimentary 
analytic work to find the critical surfaces for vanishing tension and mass gap. This could probably be 
extended, order by order in processes which change the number of particles, to derive an approximate 
scaling trajectory at intermediate lattice spacing. In practice, however, it may be just as well to 
do this numerically. Basically we are dealing with a very discrete problem in the transverse lattice 
direction and numerical methods are much more efficient. In the longitudinal directions on the 
other hand, where we take the continuum limit, we have demonstrated that semi-analytic methods 
produce a vast improvement in the removal of discretisation errors. 

An 'equal-time' Hamiltonian investigation using colour dielectric variables Q would also be 
interesting as a comparison. However, even in the large N limit, the quantum vacuum becomes 
complicated, making study of the spectrum more difficult. 

The most obvious generalization of our transverse lattice work is to 3 + 1 dimensions [Q. The 
same techniques are applicable here, although the numerical computation becomes more difficult. 
Although one must add more terms to the effective potential to a given order in the fields, there are 
many more physical observables which can be used to fix the associated coupling constants. 
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A Eguchi-Kawai Reduction 

Large- iV Eguchi-Kawai reduction at = means that the hght-cone Hamiltonian in the basis 
of Wilson loop Fock states built on |0) ( pO| ) is the same as the Hamiltonian for the corresponding 
problem where the transverse lattice is compactified on m links, m e {1,2, . . .}, id est one makes 
the identification M^i{x~) — M^i^„ii^{x^) in both basis states and Hamiltonian. In particular for 
m = 1 this means the identification 



For the basis states themselves the identification ( [77| ) is obviously a 1 — > 1 mapping since any single 
Wilson loop for a state of = is completely specified by its shape. The shape is completely 
specified by the sequence of orientations of the link modes in the colour trace and does not require 
knowledge of the absolute positions of the modes on the lattice. 

Now the original Hamiltonian (|2^) is quartic and local in the transverse direction. Therefore 
link modes on the Wilson loop interact only if they are within two lattice spacings of one another. 
Under reduction (|7^); modes which were not within two lattice spacings now lie on the same link 
and will interact. Also modes which were not on the same link, but were within two lattice spacings 
of one another, will have additional forms of interaction under reduction. The problem is to show 
that the extra interactions which arise from the reduction ( [77| ) are suppressed by 1/N. This is a 
finite task since the Hamiltonian, being quartic, involves no more than four link modes in initial 
plus final states. The full proof is nevertheless quite tedious, so we illustrate the steps with a couple 
of examples. 

Our first example is a 2 — )■ 2 amplitude which occurs on neighbouring link modes before reduction. 
It is the combination 



M^i{x-) = M{x-) for aU x^ . 



(77) 



Tr 



x^)a\i{(l' , + a)a+iik', x^ + a)a+i(fc, x^) 



q + q' ^k + k\ 



(78) 



which comes from the operator coupling to A2 in the Hamiltonian 




(79) 



On the successive link modes 




(80) 



the result is 




(81) 
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In the reduced theory the label is dropped in Eqns. (|78|), ([79|), and (pO|). The result this time is 

A2Tr{...a;i,^.(q')«li,,,(<z)---}|0> 

+^ Tr {aii(g)aii(g')} Tr {■ ■ ■ S., ■ ■ •} |0) . (82) 

The reduced amplitude picks up an extra term which is suppressed by Note that the longitu- 

dinal coordinate plays no role in the calculation. 

Our second example is a 3 ^ 1 amplitude involving widely separated links before reduction. It 
is the combination 

Tt ^a^_^^{ki,x^)a+i{ki,x^)a^i{k2,x^)a+i{k3,x^)^ , k^ = ki + k2 + ks , (83) 
which comes from the operator coupling to Ai in the Hamiltonian 

^Tt[m^iMIm^iMI} . (84) 

On the separated link modes 

Tr {• • • al,^^^{k,,x')al,^^,{k2, x') . . . al,^i„^ik3, i/) ■ ■ •} |0) (85) 

the result is zero since (^) is local to the link x^. In the reduced case however the result is non-zero 

AiTr{4;...}Tr{...a;i,„„(fc4)---}|0) (86) 

Once one takes into account the A^-normalisation of states, see (|3l|), one sees that ( p6| ) is 
suppressed. The generalisation of the proof to two or more transverse dimensions is straightforward. 

Eguchi-Kawai reduction was proved for Wilson's lattice gauge theory using the equivalence of the 
Dyson-Schwinger loop equations [ p^ . This equivalence required that, in the reduced formulation, 
extra terms proportional to the expectation value of a pair of disconnected Wilson loops with opposite 
non-zero winding number vanished. This was true under the assumption of large- factorization 
and unbroken U{1) winding number symmetry. However this U{1) is spontaneously broken at weak 
coupling (id est small lattice spacing) spoiling the equivalence under reduction. Our proof of 
equivalence in the light-cone Hamiltonian framework based on the dielectric variables in the vacuum 
|0) makes no reference to the winding sectors at all. The extra disconnected loops which give rise 
to the problems in the Dyson Schwinger approach have no analogue here because the large- light- 
cone Hamiltonian propagates loops without splitting or joining them. This is a special feature of 
the light-cone quantisation which does not hold in the equal-time quantisation. 

Does this mean we avoid the breakdown of Eguchi-Kawai reduction? The answer is yes, but 
unfortunately not in a way that that can take us to small lattice spacing. Our proof is based on the 
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dielectric variables acting on the vacuum |0). Provided this is the true vacuum, our quantisation 
is valid and we can use the reduced formulation. If the true vacuum lies elsewhere, which would 
be signalled by tachyons in the spectrum (see Fig. then even our quantisation of the unreduced 
theory about |0) is not valid. In particular, our quantisation breaks down for < 0, which puts a 
lower limit on the lattice spacing we can achieve on the scaling trajectory; see Fig. 0. At present we 
do not know how to formulate the light-cone Hamiltonian theory in terms of variables other than the 
dielectric ones and, even then, anywhere other than about the M = classical solution. Thus, even 
though the light-cone formulation avoids the usual source of breakdown of Eguchi-kawai reduction, 
we are unable to exploit this to reach small lattice spacing in the reduced theory. 

B Tamm-Dancoff Extrapolation 

Our numerical calculations of the glueball spectrum use a cutoff on the maximum number of link 
variables p in a state, id est a Tamm-Dancoff cutoff on the Hilbert space. Although a popular 
truncation scheme there is very little work in the light-cone literature on how one expects this to 
extrapolate to p — > oo. In the following, we will use the "transverse only" model of Section 4.1 to 
estimate the convergence in particle truncation in the limit of large p. 

Let us symmetrise our states |(Ti, . . . , ap) under orientation reversals 0:17;—+ crp_i+i, reflections 
T'l : (Ti — > — CTi, and cyclic permutations. The total number of distinct states, after symmetrisation, 
is bounded above by {J/2) ^^"^ bounded below by (p/2) / (^^') ■ ^^^S^ the actual number of 
distinct states should approach the lower bound because a generic many-particle state will not be 
symmetric under the above 3 operations. Thus, in the limit of large p the number of states is 

op-l 

(87) 



and the number of states for each value of 7 is 

2P+l/2-4p(7-l/2)^ ^/J^ 
9 

irp'^ 

Note that the distribution becomes increasingly peaked at 7 = 1/2 as p is increased. The normali- 
sation for a generic state is 

l = (ai,...,ap|ai,...,cTj,) = C2^,...,,^. (89) 

That is, we expect only one nonzero contraction for a generic many-particle state. Likewise, the 
number of p — 2 particle states that couple to a given p particle state is typically p/2. Assuming 
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only one nonzero contraction, the associated matrix element is, on average, 



((71, . . . , tTpl P^,P^ Icti, . . . , (Tp_2) 



p(l + 2m/3 + l/2mi + 4/^) _ 
2™r(l + 2/3)'r(3/2 + 3/3) 

(90) 



note that these two states must have approximately the same value of 7. Using perturbation theory, 
the amplitude for a generic p-particle basis state in a low energy eigenstate is 

Xcti, ...,(Tp\ P^P'' \ai, . . .,ap-2)^^^^ 



dp oc 



(91) 



2{<Ji,...,<jp\P^Pt^\ai,...,<jp) 
times an undetermined power of p (which we will ignore). Again using perturbation theory, one 
finds for a p-particle truncation that a low energy eigenvalue converges as 



2P- 



((Tl, . . . ,f7p| Pf.P'^ \ai, . . .,<Jp+2) 



(92) 



/27rp3/2 (cti, . . . , C7p+2\ Pfj.P'^ ki, . . . , crp+2) 

where we have estimated the number of states from Eqn. (|87|). Since we know from Eqn. (88) that 
the distribution of states is strongly peaked at 7 = 1/2, we set 7 to this value. (This is most valid 
when 2Ai ~ A2.) Thus, the convergence in particle truncation is exponentially fast 



Mi 



(93) 



where 



= In 



2[Ai + f]r(/3 + i)'r(2/? + i)' 



r(3/3 + i) M^^ivd + 4/3) + Ai + r(/3 + ly + ^•^'Jr(V;+l) 



27ra/r(2/3+ 1) 



(94) 



This is the formula for convergence in particle truncation that we use in our numerical work. As 
illustrated in Fig. it agrees well with numerical results for the "transverse only" model. Conse- 
quently, our biggest assumption is that it can be applied to the full theory as well. 

C Improved matrix elements for DLCQ 

In this section, we will describe a novel technique for improving the numerical convergence of calcu- 
lations which employ Discretised Light Cone Quantisation (DLCQ) |Q to solve the associated set of 
coupled integral equations. An important advantage of DLCQ in a many particle calculation is that 
the resulting Hamiltonian matrix is quite sparse and it is essential that we preserve this property 
when making any improvements. 

The integer- valued "harmonic resolution" K denotes the coarseness of the longitudinal momen- 
tum discretisation. Naively, one would expect that the leading finite K error associated with a 
DLCQ calculation would go as 1/K^\ this is simply the error associated with approximating an 
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Figure 11: Convergence of eigenvalues InA = ln|Mp — vs particle truncation p for the 

"transverse only" model. Here, we have chosen fx^ = 0, Ai = —Q.lg'^N , and A2 = —0.2g'^N. The 
points are from a numerical calculation using the basis (^3|). Also shown is a numerical fit to the 
function —0.698 + cp + 1.597 Inp where c = —1.145 is given by Eqn. (p^). 

integral using the trapezoidal rule. However, our Hamiltonian contains several singularities which 
make the convergence substantially slower. First, there are (integrable) endpoint singularities which 
go as 1/Vk+ and produce a leading error of 1/y/K. Second, there is the (non-integrable) singularity 
associated with the mass term and the instantaneous interaction 1/ (fc+ — p''")^ which gives a leading 
error of 1 /K"^^ js^ . Finally, the error for producing a state with more particles goes as 1/ \/K. 

Our method is the following: when calculating the matrix elements of the Hamiltonian in a many- 
particle calculation, we will use ordinary DLCQ commutation relations to calculate contractions 
associated with spectator particles. To calculate the interaction itself, we take the DLCQ momenta 
of the interacting particles, project onto a smooth wavefunction basis (continuous longitudinal mo- 
mentum), and calculate matrix elements of the interaction in this smooth wavefunction basis. Our 
modification of DLCQ essentially eliminates the and errors mentioned above. The 

remaining errors include a (small) 1/K error from the instantaneous interaction 1/ (fc''" — p"*")^, 1/ 
errors from DLCQ itself, and a 1/K error from production of particles at small momentum. 

We will discuss the (2 particles) — > (2 particles) and the (1 particle) — > (3 particles) interactions 
separately since we use somewhat different techniques for each. 

C.l (2 particles) (2 particles) interactions 

Consider a typical 4-point interaction 
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with longitudinal momenta fe; G {1/2, 3/2, . . .} and K — ki + k2 — kj, + k^. We start by calculating 
the first K cigcnstatcs of the two-particle bound state equation 



M'^{x) = I dyH{x,y)^{y) 



f\ ^(?/) r {x + y)i2-x-y) ^ \ 

x{l - x) 27ra Jo ^ Jx(l ~ x)y(l - y) \ (x - yf ) 



(96) 



^Jx{l - x)y{l - y) 

where a principal value prescription is assumed for the integral. Here, we have chosen 2Ai-t-A2-l-2A3 = 
—g^N so that, at /i^ — 0, the lowest wavefunction is constant with eigenvalue zero. Then we write the 
incoming and outgoing DLCQ states as linear combinations of these K eigenfunctions and evaluate 
matrix elements of the interaction V{x, y) in the eigenfunction basis. The result is expressed as a 
Taylor expansion in /3. Details are described in the following. 

We introduce a basis of even and odd polynomial wavcfunctions (labeled by subscripts s and s'): 

Hx)^{x{l~x)f+^-'\^ , (97) 

s e {1,2,..., S}, S ^ K/2, and /3 was defined in Eqn. (|4l|). We define, for both even and odd 
wavefuntions, matrix elements of the interaction term Vg^s' — ^ dx dy (j)s{x)V{x,y) (j)s' {y), the in- 
ner product matrix Eg.s' — J dx dy (f>s{x) (f>s' (y), and the two-particle Hamiltonian (|9^) Hg.s' — 
J dxdy (j)s{x) H{x,y) (j)si{y). We can Taylor expand these matrices in powers of f3: Vg^s' = Vg^Ji + 
PK^S' + • • •> Eg^g, = eI"1 + pE^2, 4- • • •, and Hg^g. = H^}, + /3H^^]l + • • •. Using the formulas in 
Refs. 11, |o| along with the definition of the Beta function B(a,/3) = da;x"-i(l - x)'^-\ vj"-*, 

(n) (n) 

E^ J, , and H], J, can be determined analytically. In addition, we define a K x S matrix Ti^g which 
maps the wavefunction basis onto the DLCQ basis (the subscript i e {1,2, . . . K} denotes a two 
particle DLCQ basis state with momenta ki = i — 1/2, k2 = K — ki) 

T.^s = T}^ + PT}^ + = (^)'"" { ^ ■ (98) 

However, the polynomial basis that we introduced has a fatal flaw: the inner product matrix Eg^gi is 
very poorly conditioned ]4l| ]. Instead, we find a new, well-conditioned, basis (denoted by subscripts 
v, p, and a) in which i?*-"-* is diagonal. Thus, for both even and odd wavcfunctions, we define the 
Cholesky decomposition |^|^ -E'"' = LL^ and calculate the matrix products 

^i.-; = iL-').^gVS{L-lg,,, (99) 

- (100) 

- {L-')^,H^l{L--)g,,, (101) 
= i;^:'(L-^)_. (102) 

The matrix elements of the lower triangular matrix are related to the coefficients of the Legendre polynomials 
defined on the interval [0,1]. 
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Because the polynomial basis is poorly conditioned, all of the above calculations, including the 



matrix products (99) - (102), arc numerically unstable and must be performed either exactly or 
using high precision arithmetic. 

Working in the well-conditioned basis, we solve generalised eigenvalue problem Hxt — htExt with 
eigenvalues ht = h^^'' + fih^^ + ■ • • and eigenvectors Xt — xf'' +PXt''' + • ■ ■■ First, we solve the leading 
order eigenvalue problem H^'^^'x^'^ = ^^t^^X^t^ (recall i?^^p is the identity). Nonleading powers in (3 are 
calculated using Rayleigh-Schrodinger perturbation theory extended to the generalised eigenvalue 
problem. Thus, if we write 

Xt = xf'+PY.^lix^' + --- (103) 



then 



h^P = xf[H^')~hT^Ei^))xT' (104) 



= ^ ,(0) ,(0) ^ ' (105) 

hi ' - hu' 

cS = -\x'fEi^\f^ (106) 



Using the lowest K eigenvectors, Xt: we define an eigenstate basis (denoted by subscripts t and t'), 

T^^t = T,,, (xt), . (107) 

The matrix Ti^t maps the eigenfunction basis onto the DLCQ basis and Ti_tTt.j is nearly equal to 
the identity. In order that the map between the eigenfunction basis and the DLCQ basis preserves 
eigenvalues, we find an orthonormal matrix that is nearly equal to T^^f . Thus, we perform a standard 
QR-factorisation T = QR^ where Q = Q(") + /3Q(i) + • • • is orthogonal and = + 



pR^^^ + • • • is upper triangular. To obtain non-leading orders in /3, we use the relation 

q(o)^q(«) + rM'^rW-'^ ^ q(")^ ||r(") - R^°^'^ (108) 

and the orthonormality condition 

Q(«)'^g(o) Q(o)^g(") g^^^^Q^""") . (109) 



Eqn. ( |l09|) deter mines the diagonal matrix elements of Q^^^ Q^"^ We can use ( |l08|) to determine 



the strictly lower triangle of Q'"' Q^") and then we use (|l09| ) to determine the upper triangle of 

-'^-'^That is, Gram-Schmidt orthogonalisation. 
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Figure 12: Convergence of the (a) ground and (b) second excited states vs two particle trun- 

cation, Ai = A2 = A3 = for ji^ = 0. The grey points are from standard DLCQ and the black 
points are from using the improved matrix elements. In (b), the eigenvalue for the improved matrix 
elements is w 14.043 + 0.012/X + 0.0085/^2. 




Figure 13: Convergence of the (a) ground and (b) second excited states vs 1/K, two particle trun- 
cation, Ai = A2 = A3 = for /i^ — g^N/i/la) or (3 — 1/4. The grey points are from standard DLCQ, 
the black line is with the improved matrix elements expanded to order (3^ , and the dashed line is 
with the improved matrix elements expanded to order (3^. 



Q{0) Q{n)_ Finally, we use ( [l08| ) to determine E^o) . From this, we can determine Q^"-' and 



We discard and use Q to rotate to the DLCQ basis, 



14,, =Qm (x*)Ik,p (xt')pQ? 



(110) 



Vij can be calculated order-by-order in (3. To see how well the improved matrix elements work, we 
can look at the eigenvalues vs K for a two particle truncation. The convergence with K for /x^ = 



and /i = g N/ (2a) is shown in Figs. |12| and |l 
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Figure 14: Region Rk^MM ^o^' = 3/2, _fir = 5/2, and = 11/2. The dots represent the points 
{ki/ K, K) and Rk^MM t'^^ surrounding polygonal region. 

C.2 (1 particle) ^ (3 particles) interactions 

Consider a typical 4-point interaction 



(111) 



where ki £ {1/2, 3/2, . . .} and K = ki + k2 + k^. We introduce a two-dimensional space of continu- 
ous longitudinal momentum fractions; if the three outgoing particles have longitudinal momentum 
fractions x, y, and 1 — x — y, the corresponding point in this space is labeled by (a;, y). We define a 
polygonal region RkiMM surrounding (ki/K, k2/K) as shown in Fig. |lj. Next we define a basis of 
wavefunctions 



- / Ck^MM (2^2/(1 -X- y)f , (x, y) € RkiMM 







otherwise 



(112) 



(113) 



where the normalisation constant C'ki,k2.ka is given by 

1=1 dx I dy(t)ki,k2-MiX'yf ■ 
Jo Jo 

We calculate improved matrix elements associated with interaction V{x,y, 1 — x — y) using the 
expression 

^fei,fe2,fe3 = / dx dycjik^M-kA^^yWix^y^'^'' X ~y) ■ (ii4) 

Jo Jo 

yk1.k2.k3 is calculated numerically order-by-order in (3. 
C.3 Implementation 

In our numerical work, we expand the interaction terms to order /S^, introducing an order error. 
Thus, for large convergence with K is somewhat degraded although the eigenvalues do still 
converge correctly. From Fig. ^ we see that this is not a significant problem if fi^ < O.bg^N/a. 
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Even though the mass term is a two-point interaction, we calculate its matrix elements as if it 
were a four-point interaction. Note that the expansion of the mass term in powers of /9, has an order 
1//3 contribution. 

In the (2 particles) (2 particles) case, we use = 40 states in our polynomial wavefunction 
basis; calculations in this basis were performed with Mathematica using exact arithmetic (except for 
Ti^s where high-precision floating point arithmetic was used). For, the (1 particle) — )■ (3 particles) 
calculations, we used standard numerical techniques to perform the associated integrations. 

In practice, we calculate matrix elements of the various interaction terms beforehand and store 
them in a file. In a many-particle calculation, the computer program simply looks up the appropriate 
matrix element in a table E2| . 
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